Controlling the composition of a confined fluid by an electric field 
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Starting from a generic model of a pore/bulk mixture equilibrium, we propose a 
novel method for modulating the composition of the confined fluid without having 
to modify the bulk state. To achieve this, two basic mechanisms - sensitivity of the 
pore filling to the bulk thermodynamic state and electric field effect - are combined. 
We show by Monte Carlo simulation that the composition can be controlled both 
in a continuous and in a jumpwise way. Near the bulk demixing instability, we 
demonstrate a field induced population inversion in the pore. The conditions for 
the realization of this method should be best met with colloids, but being based 
on robust and generic mechanisms, it should also be applicable to some molecular 
fluids. 

I. INTRODUCTION 

The question of how the adsorption of foreign particles affects the properties of materials 
and the means to control this is of central importance in domains ranging from separation 
processes to nanotechnology. This motivates the continuing investigation on the factors 
determining the adsorption process and the search of conditions most favorable for its control 
(composition of the adsorbing fluid, adsorption geometry etc.) The purpose of this 

communication is to propose a novel method that allows this, on the basis of Monte Carlo 
simulations of a fluid-slit pore equilibrium. To avoid having to consider specific interactions 
as in molecular adsorption, we choose to point out the basic mechanisms on a simple model 
with only hard-sphere and dipolar interactions. The situation closest to this model is then 
the adsorption of macroparticles. Another reason is the recent development of studies of 
colloidal adsorption Indeed, while numerous studies exist on molecular adsorption (see 
for e.g. 3, 4] and refs for more recent work), one practical advantage of using colloids is 
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the possibility to tune their effective interaction (eg. by adding polymeric depletants) and 
their coupling with external fields, possibly in confined geometry 6j. The actual behavior 
may, however, be complicated by the interplay of different effects (see for example the role 
of static and hydrodynamic forces in microfluidics devices ij]). We thus propose here a 
method that allows a fine control of the composition of the adsorbed fluid, while remaining 
very simple. 



II. METHOD 



We start with the simplest confinement geometry: an open slit pore with parallel walls 
in equilibrium with a bulk fluid. It has been used in several theoretical studies to determine 



B, y, y, q. 
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the parameters affecting the behavior of the confined fluid, (see for example 
12 . 113I . [ijj] and references therein). Since we seek a method in which the p 
the interactions (between the particles and the particles and the confining medium) as well 
as the bulk thermodynamic state are fixed, one alternative is the coupling with an external 
field. This should always be possible since besides particles having a permanent dipole such 
as magnetic colloids, colloidal particles are always polarizable to some extent. We thus took 
a uniform electric field E = E u z normal to the walls. As in [ij]] we considered a mixture 
in which one species bears a dipole moment /x, taken permanent for simplicity. The field is 
then not applied in the bulk. We thus have pure hard-spheres (species 1) and dipolar hard- 
spheres (species 2), possibly with a non-additive diameter o"i 2 in the potential ( see a ls° 
the discussion of figure 5). Both species have a hard-sphere interactions with the walls. This 
makes the model more appropriate ll| to a mixture of hard-sphere-like colloidal particles, 
than to a molecular mixture (see however the final remarks). The effect of an external field 



Sd temperature) on the filling of a cylindrical pore was also studied in 16] (see also and 
I for a slit pore), but not from a bulk mixture. Previous studies considered the role of 
n 

;he pressure in one-component fluids (eg. |8(), or the total density and the mole fractions 
9), Q] in bulk mixtures but without field. As shown below, the combination of both 
will play here a crucial role. An inhomogeneous multicomponent mixture with anisotropic 
interactions being difficult to study by analytical methods (we are aware of one study by 
density functional theory of the adsorption from a mixture of polar molecules [3]), we used 
Monte Carlo simulation (see also, for e.g. [l5, 16, 19, 20, 21]). We already pointed out how 
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the structure can be modulated by the combination of various interactions 10|, L3J| arid by 



the action of the field 
could be modulated in 



15 1 . However, only the density profile of the particles through the pore 
151 ] since the total number of particles was kept fixed (simulations 



in the canonical ensemble). An important difference here is that we consider an open pore 
which exchanges particles with a reservoir. One may then achieve much stronger variations 
of the density of each component in the pore. The physical pore is assumed large enough that 
the interfacial region in which it is in contact with the reservoir plays a negligible role 22]. 
For this reason we will refer to the fluid in the reservoir far from this region as the "bulk". 
The pore/fluid equilibrium is determined by the equality of the chemical potentials pi and 
P2 of both species in the bulk and in the pore. But since the practical control variables are 
the total density pb and mole fraction X2 of the dipolar species, in the bulk, the latter is 
studied in the canonical ensemble. By considering only homogeneous states or metastable 
states very close to the coexistence boundary, pi and p 2 in the bulk are determined with 
sufficient accuracy from Widom's insertion method (see for eg 23J] for this point). p\ and 
fi2 are then used to study the fluid in the pore in the grand-canonical ensemble. We can 
then compute the average density of each species in the pore as a function of pb and X2- 
Hereafter, reduced variables E* = E(a 3 /kT) 1 / 2 and p* = p/(kTa 3 ) 1 ^ 2 will be used. The 
reduced density in the pore is p = Na 3 /V, with iV the average number of particles for a 
lateral surface S with periodic conditions in the x and y directions and V = S(H — a) the 
accessible volume in the pore. We took a pore width H = 3cr. In the bulk, iV = N\ + N2 is 
fixed. 



III. RESULTS 



We show in figure 1 the first basic mechanism: a field induced filling of the pore by a 
one-component dipolar fluid. At increasing field strength E*, the pore is progressively filled 
by the dipoles with a rate that depends on pb . For the value of the E* and p* used here, 
the explanation seems that the field-dipole interaction energy — J2i Vi-E offsets the entropy 
loss due to their orientation in the direction of the field. This is more visible at low pb in 
which case the slope is nearly constant beyond E* = 8. For a particle diameter of lpm and 
T = 300-fT, for example, this corresponds to E = 49 10~ 3 V/ pm and p = 2 10 5 D. Thanks 
to the scaling factor a~ 3 l 2 in its definition, the same value of p* is also appropriate for the 
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dipolar interaction between molecular species. 




10 12 14 



FIG. 1: Effect of the applied field on the filling of the pore by a one-component dipolar 

fluid. 

p is the total density in the pore and E* the field strength in reduced units. The bulk 
density is pb = 0.51 (filled squares) and pi, = 0.0102 (empty squares). The lines are a guide 

to the eye. 



The dipole moment being then of the order of one Debye, a field strength of the order 
10 V /nm is needed to obtain the same reduced energy —p*E*. Such field strengths are not 
unusual for confinement at the molecular scale (for water in nanopores see for example [l6] 
and [l7], in particular figure 1 in the last one). Note that the equilibrium state in presence of 
the field may not always be the filled one at other parameters or if more complex interactions 
are considered (see ref lj| for example), due to the competition between energetic and 
entropic contributions. 

Figure 2 illustrates the second mechanism: the relative population of a pore in equilibrium 
with a mixture having a natural tendency to demix. The simplest one is the mixture of non- 



additive hard spheres 



24l | in which the cross diameter is = 1/2 (c7j + o~j)(l + 5). Previous 



studies have shown that when the pore is in equilibrium with a mixture in which 

one species is in minority (say X2 = 0.02 for a non additivity parameter 5 = 0.2 ) a population 
inversion occurs in the pore when the total bulk density is varied. This occurs here for pb 
between 0.55 and 0.56 for a pure non-additive HS mixture and between 0.54 and 0.55 for 
the hard- sphere dipole mixture. 
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(a) Mixture of symmetric non-additive (b) Mixture of hard-spheres and dipolar 

hard-spheres. Filled circles: adsorption; Open hard-spheres with S = 0.2, fi* = 1 and x<i = 0.02 ; 

circles : dcsorption; the bulk concentration of the filled circles: adsorption of the dipolar hard-spheres, 
adsorbing species is = 0.02. The non-additivity 
parameter is S = 0.2 



FIG. 2: Population inversion in a pore in equilibrium with a bulk mixture close to 

demixing. 



Having the basic ingredients, we may now combine them to produce the desired effect: 
by choosing the composition of the bulk fluid so as to be close to the population inversion in 
the pore, we anticipate that the closer we are from the threshold density the weaker will be 
the external field E\ r required to trigger it. This is shown in figure 3. In the most favorable 
case shown the actual value of E tr is about 3 10~ 3 V/ \im. 

A small variation of the applied field produces the inversion: the dipolar particles are 
selectively absorbed by a weak change about E tr , the converse being possible, perhaps 
with some hysteresis 13|. Near the adorption jump (p& = 0.53, E* = 0.5 for p* = 1) 
a slight change in temperature (5T = 25°C in figure 4) produces a detectable change in 
adsorption. This observation may be important for some applications (since the reduced 
variables combine E and fi with T, while E = in the bulk, one has to rescale E* by the 
factor T 1 / 2 to compare two temperatures at a given values of E). 

Finally, in order to check the effect of a different interaction w 12 between species 1 and 2, 
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FIG. 3: The curves show the dipoles mole fraction in the pore as a function of the reduced 
applied field strength. The curves are for a bulk mole fraction x 2 = 0.02 and bulk densities 
(from right to left) pb = 0.51, 0.53, 0.54. The inset shows the corresponding dipoles and 
hard-spheres density in the pore for pb = 0.53. 
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FIG. 4: Adsorption-field strength curves at different temperatures. 
Dipoles mean density in the pore as a function of the applied field strength at T = 285K, 
300-fT and 325K from left to right. The bulk mole fraction is x 2 = 0.02 and bulk densities 
are p b = 0.53. Here E*(T) = E*(300)(T/300) 1 / 2 . 
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we also show in figure 5 the result for a Yukawa repulsive potential. 



0.8 



0.6 



0.4 ■ 



0.2 - 



1 1 


i i 


1 1 1 


i 

„ ■ 


h 


/ f 

I / 

' / 






/ 
/ 

/ 
/ 

/ 
/ 

/ 


/ 

/ 

/ 

/ 

/ 






/ / 

// 


1 1 


i i i 





0123456789 

# 

E 



FIG. 5: Effect of the applied field on the filling of the pore by different model fluids. 
p is the density of the dipolar species, in the pore. Empty squares : one-component dipolar 
fluid (for pb = 0.0102 as in figure 1); open circles : additive mixture of hard spheres and 
dipoles (with x 2 = 0.02 and pb = 0.51 in the bulk). Diamonds : same for a Yukawa 

repulsion between the dipoles and the hard-spheres; filled circles : same for the 
non-additive hard-spheres - dipolar hard-spheres mixture. The range of the Yukawa 
potential (with e* = 8) gives the same contribution to the second virial coefficient as the 
non-additive hard-spheres potential with 8 = 0.2. Filled squares : one-component dipolar 

fluid with pb = 0.51 (as in figure 1). 



We observe that the phenomenon is quite general. The requirement for observing a 
sensitive field effect is that the self coordination should be more favored in the mixture. The 
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ity can be favored by suitable chemical composition of the particles surface 
27]] . Weak specific interaction with the pore walls can be achieved similarly. 
As our main goal was to demonstrate the phenomenon of a field activated adsorption from 
an unstable mixture, we used a simple simulation strategy. Accordingly, we did not conduct 
a detailed study of the behavior of the confined fluid. For instance, phase equilibrium in 
the pore may take place before the spontaneous condensation (adsorption jump) predicted 
in the grand canonical simulation 25 j - According to Duda al. jl^. the inversion line for 
non-additive hard spheres is close to the bulk fluid coexistence line but the two phenomena 
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are different. We actually observed that the inversion corresponds to a bulk fluid close to 
demixing or slightly in the two-phase region. Regardless of this, the essential point is that 
the density in the liquid-like phase should be close to the value after the adsorption jump, as 

The precise relation between these observations and other 



m one-component systems 



phenomena such as capillary condensation, wetting, hysteresis, etc. (see e.g. 25 
will be discussed in future work. 
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IV. CONCLUSION 



In conclusion, these results show that the combination of two generic mechanisms allows 
a quite sensitive control of the pore filling. Although this method has been demonstrated 
for particles that are closest to the optimum conditions (ie hard-sphere-like colloids), none 
of them is exclusive and since the basic mechanisms (demixing instability and coupling with 
an external field) are quite generic, this prediction should concern a broader class of systems 
(including molecular ones). In order to benefit from the field effect, one species should be 
either polar (e.g. ferrocolloid in magnetic fields) or much more polarizable than the other (the 
results given here being relative to permanent dipoles). The solution should also not contain 
free charges to avoid particle motion due to the action of the field (electrohydrodynamic 
flows), not considered in this simple model. Polar molecules being on the other hand rather 
common, one should consider in this case also the role of specific interactions. We believe 
that further experimental studies and simulations on this method are worthwhile given the 
diversity of possible applications of this field controlled composition of the confined fluid 
and hence flexible control of the physical properties that depend on the composition of the 
confined fluid. Just as an example, one may consider to modulate in this way the dielectric 
response of the confined fluid for optical applications. 
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